00 

o 

O 
(N 



Tj ; 

o ■ 

Oh. 



00 

o 

oo 
O 



X 



A Lorentz and gauge invariant measure of laser intensity 
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Focussing on null fields as simple models of laser beams we discuss the classical relativistic motion 
of charges in strong electromagnetic fields. We suggest a universal, Lorentz and gauge invariant 
measure of laser intensity and explicitly calculate and interpret it for crossed field, plane wave and 
vortex models. 
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I. INTRODUCTION 



Since the advent of chirped pulse amplification [l| op- 
tical lasers have reached unprecedented intensities, the 
present state of the art being about 10^^ W/cm^. The 
associated electric fields of approximately 10^^ V/m have 
renewed interest in the idea of particle acceleration using 
laser fields. Using plasma wake fields impressive progress 
has been made during the last few years [4]. In vacuum, 
however, achieving laser acceleration is nontrivial as a 
laser field is oscillatory; the periodic sign change seems to 
render acceleration of charges unfeasible. Indeed there is 
a result named the Lawson- Woodward theorem [3| which 
states that a plane wave cannot accelerate charges (see Q 
for a thorough discussion) . This is somewhat academic to 
the extent that a laser beam is only crudely modelled by 
a plane wave, as there is no ponderomotive force. Nev- 
ertheless, adopting plane waves as a first approximation, 
it suggests that a laser may not be a very efficient parti- 
cle accelerator in vacuum - unless its oscillatory and/or 
plane wave characteristics can be overcome or modified. 
Part of this paper will be to analytically assess some op- 
tions in this direction. 

The most important parameter characterising the 
'strength' of a high power laser is the 'dimensionless laser 
amplitude' which, in a particle physics context, may be 
written [5[ 



eSri 



ao 



LorrieC 



(1) 



with e being the elementary charge, E'rms = {E"^)^/"^ the 
rms electric field, u the laser frequency and me the elec- 
tron mass. The analogue of 1/ao in atomic physics is the 
Keldysh parameter, see @ and references therein. Em- 
ploying the laser wavelength Xi^ = c/uj we may write 



eEr, 



Oo 



(2) 



which describes ao as the ratio of two energies; the aver- 
age energy gain of the electron in the laser moving over a 



laser wavelength, divided by its rest energy. The presence 
of the velocity of light, c, signals the relativistic nature 
of the strength parameter, with ao > 1 describing the 
regime where electrons become relativistic. 

Note that ao is a purely classical quantity as no fac- 
tors of h are present. This should be compared with 
the situation when the denominator of ([2|) is written 
in terms of Sauter's critical field [7|, more commonly 
known as the Schwinger limit [8|. This is defined by 
rrieC^ — eEcrit'^Ci where Ac = h/rrieC is the Compton 
wavelength. At the critical field strength it becomes en- 
ergetically favourable for the vacuum to 'break down' by 
emitting electron positron pairs which will move apart 
such as to decrease the external field until it becomes 
subcritical again. This is clearly a quantum effect. Nu- 
merically, one finds Scrit = 1.3 x 10^^ V/m, which is 
still five orders of magnitude above what can presently 
be achieved. Inserting this definition into ([2|) we find 
flo = {^L/^c){Erms/Ecrit)- The smallucss of the second 
factor is compensated by the large ratio of laser to Comp- 
ton wavelength which is of the order of 10^. That is why 
the high-intensity regime, which is the main focus of this 
paper, indeed corresponds to ao > 1, the current upper 
limit being about 10^. 

Looking at (H]) we note that ao appears to be a frame 
dependent quantity as the laser wavelength (and, more 
subtly, the time average involved in defining -Brms) do not 
represent Lorentz invariant concepts. This shortcoming 
is claimed to be overcome |5[ by employing the gauge 
4-potential A^^ to define ao via 



e^{A^A^^) 



(3) 
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This new definition, however, seems like a cure worse 
than the disease as (for a properly chosen average) it es- 
tablishes Lorentz invariance by sacrificing gauge invari- 
ancc - ([3]) is not equal to ([1]) in all gauges. This is clear 
from a particle physics perspective as A^^A^^ precisely cor- 
responds to a photon mass term (say in a Lagrangian) 
which certainly violates gauge invariance. 

The main topic of this paper will hence be to provide 
a universal definition of ao that meets both the require- 
ments of Lorentz and gauge invariance. Adopting com- 
mon particle physics practice we will set h = c = 1 in the 
remainder of this paper. 



II. 



MOTION IN PLANE WAVE FIELDS first integral of the equation of motion (g]) is 



Since we know from ([T]) that ao is a classical quan- 
tity, it will suffice to consider the classical motion of a 
charged particle (chosen to be an electron throughout), 
as governed by the Lorentz equation 



p^{r)^—F,,{x{r))p-^{r), 



(4) 



maintaining both gauge invariance and Lorentz covari- 
ance. Throughout, p^ = niei'^ is the 4-momentum of the 
electron and dots represent derivatives with respect to 
proper time r. The right-hand side of ([4]) is the Lorentz 
force acting on the particle. Gauge invariance is mani- 
fest as only the gauge invariant field strength tensor Ff^^^ 
appears^. 

Our first task is to solve (j4]) in terms of F^i/ rather than 
the potential Afj, (recaU F^i, = dfj,Ai, — dj^Ai^) which was 
the previously adopted approach [l^, [ll| - see, however, 
[12l | . For an arbitrary field configuration it will generally 
not be possible to find an analytic solution of ^ which is 
a nonlinear differential equation depending arbitrarily on 
x(t) via -F)i^. On the other hand, drastic simplifications 
arise if F^y{x) is a plane wave. In this case the field 
strength tensor may be written [8| 



F^'^ix) ^F,{k- x) /r + F2(fc • x) /r 



(5) 



with the amplitudes Fj depending on the Lorentz invari- 
ant phase k ■ x, where k'^ = ujn'^ is the 4-momentum of 
the wave. The constant tensors /'"' are parameterised in 
terms of dimensionless 4-vectors n and e 



P^{t)^cxp[G{t)]^^p'^{Q) 



(7) 



To evaluate the exponential in ([7]) we have to determine 
all powers of F^^ (or G^^). For F^^ one finds the rather 



simple result. 



Ci2 = p crpi 



FjFj n^r 



(8) 
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Multiplying with F^j^ once more, transversality implies 
that its cube vanishes, F^^^ = 0. Fields with t his p rop- 
erty are called null fields in general relativity [ij|. In 
Taub's classification of constant fields this case is de- 
noted 'parabolic' [IJl i-e. electric and magnetic fields are 
orthogonal and of equal magnitude such that the basic 
bilinear Lorentz and gauge invariants vanish. 

-iF^.F^'^ = i(F'-S') = 0, (9) 

-^F^.Ff^" ^E-B^O. (10) 

Constant fields of this type are also referred to as crossed 
fields, e.g. in [11|, and constitute the simplest model of 
a laser beam applicable when all relevant length scales 
are small compared to Xl- As (l9|) and (fTO|) are invariant 
statements, identifying fields as null is a Lorentz invariant 
characterisation. It is easy to check that these identities 
are equally valid for the plane waves Q . We note in this 
context that the energy momentum tensor for null fields 
takes on a particularly simple form and coincides with 
([5]) as the trace term vanishes, T^i, = F^^^. 

As a result of this discussion the exponential series ([7]) 
is truncated to second order, giving the simple expression 



/^ 1/ 



n^e 



n e- 



(6) 



The lightlike vector n^ may be chosen 'orthogonal' to 
the spacelike polarisation vectors e^, i.e. n ■ tj = 0. For 
circular polarisation we have ei • £2 = and for linear 
polarisation we simply set one of the Fj to zero. In ei- 
ther case we clearly have n^F'^'^ ~ 0, which expresses 
the transversality of the wave. Hence, dotting n (or k) 
into the equation of motion we find that the momentum 
component k-p is a constant of motion which, upon inte- 
grating, induces a proportionality between proper time r 
and 'light-cone time' n ■ x (assuming n ■ x{0) = 0). For 
the phase k ■ x in ^ this implies k ■ x = fir. Thus, quite 
remarkably, we can trade the x dependence of F^"^ for a 
dependence solely on proper time r whereupon the equa- 
tion of motion becomes linear. Consequently it can be 
solved analytically by exponentiation. As the time de- 
pendence resides in the scalar prefactors Fj, the tensor 
F^^ commutes with itself at different times and time- 
ordering of the exponential is unnecessary. Denoting the 
integral of eF^i^/TTig by G^^ such that Gp,y(0) ~ 0, the 



Pf.{r) = [.9p. + G^^T) + ^GI,{t)] p'-iO) , 



(11) 



with g^i, = diag(l, — 1, — 1, — 1). This constitutes the 
gauge invariant first integral of the Lorentz equation ([¥]). 
It is a trivial matter to integrate again, which yields a 
gauge invariant expression for the particle orbit. For the 
moment, (|lip is sufficient for a proper determination of 
the strength parameter aq. 

From Q, where a§ is quadratic in A^, we expect a 
gauge invariant definition to be quadratic in F^^ or G^^, . 
Furthermore it should be related to an energy density 
or intensity, which suggests using the energy momen- 
tum tensor T^j,. As the usual Lorentz invariants ^ and 
(|10p vanish we are lead to the only remaining invariant 
PijT^'^Pu given by contracting with the gauge invariant 
momentum (jlip . There is more than one way to make 
this dimensionless, but our previous formulae ^ and ([2|) 
uniquely lead to the definition 



{k-pf 



(12) 



^ For a discussion of the gauge invariance of ao using A^ see [^l • 



Here {{■■■)) denotes the Lorentz invariant proper time av- 
erage which, in this paper, we take over all proper time 
such that it yields the Fourier zero mode. 



The invariant numerator in (|12p . originally introduced 
in [iJI as {F^^p^Y, is indeed proportional to the en- 
ergy density seen by the charge, as in the charge rest 



frame, where p = (ttIe, 0), one has p^T^^'^pi, 



2rpOO 



m1{E'^ + B'^)/2. The general relativistic analogue of this 
is an important characteristic of curve d sp ace times, con- 
strained by positive energy theorems [l5[. The numera- 
tor of p^ is then the average energy density seen by the 
charge during its entire exposure to the laser field. The 
second invariant k-p = mgO, in the denominator is the 
constant of motion proportional to the laser frequency 
seen by the particle. Altogether, oq is a ratio of two 
energies as described below ([2]). 

Taking p as the momentum of a probe photon ao also 
determines the amount of vacuum birefringence in ultra- 
strong laser fields due to an induced effective metric, 
h^u = giiu - nTf^i, with K = K(ao) [ll]. 

We conclude this section with three explicit examples: 
linearly and circularly polarised plane waves, and crossed 
fields. In each case our definition p^ recovers ([T|), and 
the gauge non-invariant expression ([3]) is obtained from 
(fT^ by substituting for F^" in terms of the gauge poten- 
tial A^ in light-cone gauge, n ■ A ~ [I4I. We will now 
describe the orbits and calculate ao for our three exam- 
ples. In all cases we choose initial conditions €j-p{0) = 0, 
i.e. the electron has zero initial transverse momentum, 
and as before n-x{0) = 0. From pT|) . we find that the 
orbits may be parameterised by three functions 61, 62, ^3. 
These are derived from the Fj describing the laser fields 
as in ^ and ([5]), and will be given below for each exam- 
ple. The form of the orbits is: 



rUe 



k-p ^ 2{k-pf 



(13) 



The first two terms represent initial conditions and the 
final two terms describe the transverse and longitudinal 
motion, respectively. We proceed to the examples. 
Linear polarisation. We have Fi = loC cos k ■ x, with 
C a constant amplitude, and F2 = 0. The orbit is 
described by 61 = C(l — cosilr), 62 = and 63 = 
C^(20r — sin2rir)/4. With transverse coordinates x, y 
and longitudinal coordinate z, we see a frequency dou- 
bling in the longitudinal orbit coefficient 63. Accord- 
ingly, and in the frame where the average drift veloc- 
ity vanishes, the motion in the x-z plane acquires the 
well known 'figure-8' shape [ll|, [l3| which is nothing 
but the Lissajous figure corresponding to frequency ra- 
tio 2:1, see Figure [T] The energy- momentum tensor is 
rj^fiv _ Q2j^iif^u ^,Qg2 ^^ ^j^j-^ ^Yic cosine averaging to 1/2 

so that al = e'^C^/2ml. 

Circular polarisation. We choose Fi as above and 
F2 = —wC sin fir. In this case 61 is as above, 62 = 
~C{i^T — sin Sir) and 63 = —2Cb2 i.e. there is no fre- 
quency doubling. In the average rest frame where all drift 
velocities vanish the trajectory is an ellipse in the x-z 
plane. The energy momentum tensor becomes constant, 
T'^" = C^k^k"^ , implying a laser amplitude ao — eC/nie. 



In both these examples our Oq defines the effective mass 
of an electron in the field [I^ via the square of the aver- 



age 4~momentum, 



mg (l-t-ag), identifying 1 



with a gamma factor squared [18|. Periodicity of the 
above fields results in periodic orbits, hence no net accel- 
eration. This is different for our final example. 

Crossed fields. These are obtained by letting fJ ^ in 
either of the cases above. We find Fi = Cut and F2 = 
implying constant field strength. Transverse and longi- 
tudinal orbit coefficients become quadratic and cubic in 
T respectively: 61 = CQ^t'^/2, 62 == 0, 63 == C^fl^r^/S. 
As the charge does not experience a full period of the 
laser field there is no cancellation between positive and 
negative field regions and we hence have net acceleration 
(subcycle acceleration). The energy momentum tensor 
and ao are the same as for circular polarisation. 



III. ELECTROMAGNETIC VORTICES 

A more sophisticated null field model of a laser beam 
is the electromagnetic vortex found by Bialynicki-Birula 
[2(| . The author states that this approximates a circu- 
larly polarised Lagucrre-Gauss beams in the vicinity of 
the beam axis/vortex line. The tensor structure is as in 
([5]) and © but the amplitude functions now depend on 
transverse coordinates as well as on the phase k ■ x. It is 
useful to introduce complex notation for the amplitudes, 
■^ = Fi + iF2, the polarisation vectors, e = ei + ie2 and 
the transverse coordinates, Z = e ■ x. 

The vortex field strength tensor 



F'"' = i .^{n^e" - n^e^) + c.c. 



(14) 



is determined by the complex 'envelope' function, ^ = 
Cuj Z exp{—ik-x) . It follows that the energy momentum 
tensor is of type ^, 

T'"' = l^pTi^i'' = C^iv^lZl^n^'n" , (15) 

but now depends explicitly on the transverse coordinate 
Z(t). Nevertheless, this form implies that, according to 
(I12p . the dimensionless intensity parameter for electro- 
magnetic vortices is 



2 e^C' 



-'0 — 2" 
mi 



{{ZZ*))^cot{{ZZ*)) 



(16) 



with cyclotron frequency tOc- To determine this elegant 
expression explicitly we obviously need to find Z by solv- 
ing the equation of motion, 

x^ = ^ujujce-x{e-x n^ — n-x e''} -|- c.c. . (17) 

We point out that this equation is nonlinear and cannoi 
be solved by exponentiation. Fortunately, Bialynicki- 
Birula has shown there is still an exact solution. He no- 
tices that longitudinal (n'') and transverse (e^) motion 
decouple. The former is easily solved because the tensor 





FIG. 1; Left: motion in a linearly polarised plane wave (elec- 
tron's average rest frame). Right: bounded motion in the 
vortex, projected into the transverse plane. 



initial position and momentum, e.g. if the particle is 
at rest initially at transverse co-ordinates xq, yo, the 
intensity parameter is 



cog xl (1 + K,)i2-K)+ yl (1 -k){2 + k) 
2 (1 - ^2) 



(21) 



This dependence on the particle's orbit is not so surpris- 
ing: the same is actually true for the plane wave, as ao 
contains e, me and an average over the particle's proper 
time. As the vortex has a non-trivial spacetime depen- 
dence, any measure of its strength should reflect this. 



structure again guarantees transversality so k ■ p is con- 
served exactly as before. This conservation renders the 
transverse equation linear. Its solution amounts, in our 
notation, to introducing a co-moving basis e,- via 

C(t) = e-'"^/2 ^(^) ^er-x. (18) 

In terms of C, the transverse equation of motion is 

C + inc~^C~ujcnc ^0, (i9) 

which is clearly linear, describing a damped harmonic os- 
cillator. Its solution is given explicitly in [201, ^^ we do 
not reproduce it here. For k = AtOc/^ < 1 it describes 
non-periodic bounded motion, see Figure[T] k > 1 yields 
acceleration which, however, is an artifact of the unphys- 
ical transverse growth of the fields. 

With the solution in hand, we return to uq. It is clear 
from (|18p that the intensity parameter becomes 



«0=^c((CC)>, 



(20) 



and it remains only to evaluate the proper time average 
over transverse coordinates. We find the resulting ao 
to depend explicitly on the charge's motion through its 



IV. DISCUSSION AND CONCLUSION 

We have given a thorough discussion of particle motion 
in null fields from a particle physics perspective. The 
fields in question are regarded as simple models of laser 
fields, for which the usual Lorentz and gauge invariants 
quadratic in electromagnetic fields vanish. It follows that 
the laser fields can only be characterised in terms of a par- 
ticle probe. Its 4-vector p then may be employed to form 
a nontrivial invariant PfjT^'^pi, which is the energy den- 
sity of the field as seen by the probe in its instantaneous 
rest frame. By making this quantity dimensionless one 
arrives at a universal definition of laser strength satisfy- 
ing the crucial requirements of both Lorentz and gauge 
invariance. For more complicated beams oq becomes de- 
pendent on transverse coordinates and hence acquires a 
profile. Its nonvanishing gradient will in turn determine 
the ponderomotive force. 
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